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Abstract 

The model recently proposed in Ref.l is used to derive linear integro- 
differential equations whose solutions provide reasonable estimates for the 
momentum distribution and condensate fraction in interacting many-boson 
system at zero temperature. An advantage of these equations is that they 
can be employed in the weak coupling regime and beyond. As an example, 
analytical treatment of the weak coupling case is given. 

PACS: 05.30.Jp, 67.90. +z 

A new way of investigating spatial particle correlations has recently been proposed for 
the many-boson system in the ground state |T[. It results in integro-differential equations 
for the pair distribution function. They accurately take into account the short-range boson 
correlations and, in the weak coupling regime, yield thermodynamics which reasonably agrees 
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with the data of the Bogoliubov approach |],|] . The structure of the derived equations is 
similar to that of the Born-Green equation [[| for the pair distribution function of classical 
liquids. This gives an optimistic view on the possibility of using methods developed for 
classical liquid in investigating the Bose one. The most important peculiarity of the new 
approach is that it does not assume a small depletion of the zero momentum state. Moreover, 
the density of the condensate particles is not at all a parameter contained by the integro- 
differential equations. So, a question arises how to calculate the condensate density within 
the model of Ref.l. In addition, it is useful to clarify how to determine the momentum 
distribution of particles in the Bose liquid with the help of the integro-differential equations 
mentioned above. Answering these questions is the aim of the present Letter. 

The phenomenon of the Bose-Einstein condensation consists in the macroscopic occu- 
pation of the zero momentum state so that for the total density of bosons we have (in the 
thermodynamic limit) 

n = n + j—^ J n(q) d 3 q , (1) 

where uq denotes the density of condensate particles, n(q) stands for the distribution of 
bosons over nonzero momenta. According to expression ([I]), no can be found with the 
known distribution n(q) at any given n. Hence, we would be able to calculate no within the 
approach of Ref.l if we managed to derive an expression for n(q) in terms directly connected 
with the spatial correlations. An attempt is natural to get this expression with the help of 
the ground-state energy E because there are two important relations linking it with both 
the pair distribution function and n(q). These relations are the consequences of the well- 
known statement which is often called the Hellmann-Feynman theorem and results in the 
following equality: 

5E = (ip | 5H | V) , (2) 

where SE and SH are infinitesimal changes of the ground-state energy and Hamiltonian, 
respectively, and ip denotes the ground-state wave function. The first relation is given by Q 
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N 2 f f 

E = E id + — J d 1 J g{r ] - f ,n)<5>{r)d 3 r . (3) 
o 

Here E^ stands for the energy of noninteracting particles {E^ = represents the most gen- 
eral case for bosons), V is the system volume and N = n V . Besides, 7 denotes the coupling 
constant, $(r) stands for the interparticle potential and g(r;j, n) is the pair distribution 
function. Note that below, for the sake of brevity, the notation g(r) is also used instead of 
g(r;j,n) so that g(r) = g(r;^f,n) . The quantities related to g(r) are handled in the same 

way. The second important relation mentioned above concerns the functional derivative of 

h 2 q 2 

E with respect to the one-particle kinetic energy T(q) = — — and is expressed as 

m) = Wf n(q) ' ,/0 - (4) 

Combining (0) and (HI) leads to the following equality: 



n(q) =A^n 2 I d 7 / $(r) »^ ' d A r . (5) 







Now, we have got the connection of n(q) with g(r). However, to employ it one needs 
to know more than simply the dependence of the pair boson distribution on r, 7 and n . 
One should also have the knowledge of the functional dependence of g(r) on T(q) . The 
integro-differential equations for g(r) proposed in Ref.l, turn out to provide all the necessary 
information. Let us consider the simplest of them, to shorten the further reasoning without 
loss of generality. This equation may be written in the following form: 

=7$(r)+n 7 [ H\r-y\) (u 2 (y) + 2u (y)) d 3 y , (6) 
m(l+u(r)) J v J 

which is convenient in reaching our aim. Here u(r) = g 1 ^ 2 ^) — 1 . Taking the functional 
derivatives of the left- and right-hand sides (l.h.s. and r.h.s.) of (j^) and, then, equating 
one to another, we are able to find an equation for 

5u(r)_ 
ql ' ~ 8T{q) ' 

and Kq(r) may be used in evaluating n(q) with the obvious relation 



8T(q) 



2(l + «(r))lf q (r). (7) 



h 2 

While calculating, one should realize that the form of the operator — V 2 is completely 

m 

h 2 

specified by the shape of T(q) . Hence, the functional derivative and operator — V 2 are not 

m 

commutative in the situation considered. Now, let us perturb the one particle kinetic energy 
replacing T(q) by T(q) +5T(q) . Working only to the first order in the perturbation, for the 
change of the l.h.s. of (|6|) we have 

2 / h 2 \ h 2 
S (lh.s.) = - 5 -— V 2 u - 5u VV (8) 

v ' 1 + u \ 2m J m(l + u) 2 v ' 

The relation 

Ti 2 i r 

~2m^ 2u ^ = (2tt)3 J T (<?) ^(<?)exp(2qr) rf3 ^' 
where u(q) stands for the Fourier transform of u(r), makes it possible to get 

6 ("1^) = " 2^ y2 {6u) + ^ky 1 6T{q) %) exp(zqr) d3q ■ (9) 

Inserting (^) into (|3|) we arrive at 

5{l.h.s.) = ,f . V 2 (6u) .f X9 6u V 2 u - 

y ' m(l + u) m(l + u) 2 

exp(iqr) 



i+w^^w" 9 ' (10) 

In its turn, the leading term produced by the perturbation in the r.h.s. of @ is given by 

5(r.h.s.) = 2n7 J $(|r-y|) (1 + u{y)) 5u(y) d 3 y . (11) 

Using relations (|i0|) and (|Tl| ) and the equality 5(1. h.s.) = 8(r.h.s.), one can derive the 
following integro-differential equation for i^ q (r): 

h 2 / 2 Kq(r) „ 2 ^ exp(iqr) 



v - tt^ v u ( r ) = 5(9) -^^ir- + 



2ml qw l+u(r) y ' ) y ' (2 



+ «7(l + «(r)) / $(| r - y |) (1 + u(y)) K n (y) d 3 y . (12) 



The boundary conditions for u(r) are given by 

lim u(r) = 0, u(Q) < oo . (13) 

r— >oo 

To know the behaviour of -Kq(r) , when r — > oo, one should keep in mind that the boundary 
conditions ([13|) are not changed by perturbing the one-particle kinetic energy. This implies 
that for any regular 5T(q) <C T(q) we have 

8u(r) = J K n (r) 5T(q) d 3 q -»■ (r -»• oo) . (14) 

Hence, one could expect that -Kq(r) — > (r — >■ oo) . However, it is not the case. To 
be convinced of this, let us investigate (^) and ([12|) in the weak coupling regime which is 
introduced with the inequalities 

u{r) < 1 , Jf q (r) < 1 . 

Keeping only the leading terms in ([]) and ( TT2| ) results in the following expressions: 
h 2 r 

— - V 2 n(r) = 7$(r) + 2^7 J ®(\ r - y \) u(y) d 3 y (15) 



and 
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2m V 2 K q (r) = %)^p + n 7 |$(|r-y|)^ q (y)^. (16) 

From (|l~5D and (|i~6"D it follows that the Fourier transforms of u(r) and Kq(r) (u(p) and K q (p), 
respectively), obey the relations 

5(P) = 4 (17) 

V 7 2 T(p) + n 7 $(p) V 7 

q ^ ^(p) + n r y^(p) ^ ^ ^ ^ 

Remark that in the weak coupling case the value u(p) at any given p depends only on the 
quantity of the one-particle kinetic energy at this very p . So, the 5-function behaviour of 
fT q (p) in ( TT8| ) is a result of, say, the "macroscopic" contribution of the T{p) to u{p) . With 
(|l~7|) and it is not difficult to get the expression for Kq(r) in the weak coupling regime: 



As it is seen, for r — > oo the quantity -Kq(r) oscillates and does not tend to zero. However, 
relation (||) is fulfilled. 

To have a reasonable idea concerning the long-range behaviour of i^ q (r), let us note 
that ([12|) is reduced to (|i~6D not only in the weak coupling regime but in the case r — > oo 
as well. This can easily be determined via inserting u(r) ~ (r — > oo) into ([12]). Then it 



is reasonable to expect that the solution of ( |I2]) is close to that of (p~6f ) at large distances. 
Hence, investigating ( ]T2]) we may adopt the following ansatz: 

K q (r) = D q (r) +C q exp(zqr) , (20) 

where the boundary conditions 

limL> q (r) = 0, D q (0) < oo (21) 

are fulfilled for the "short-range" function -D q (r) at any q. Substituting ( j20[) into (|i~2[) and 
using (|i~3"D and (^l|), we obtain 

C Q = - %) , . (22) 

q (2tt)3 T (g) + n 7 $(g) 

So, the "long-range" part in (|20| ) is very similar to the weak coupling expression for Kq(r) . 
However, u(q) is now the solution of the exact equation (|[). 

Integro-differential equations (^) and ([12]) provide the detailed information concerning 
the distribution n(q) . However, when evaluating only the condensate fraction n^/n, we do 
not exactly need the function -Kq(r) . It is sufficient to work with the quantity 

K(r) = J K n (r) d 3 q . (23) 

Indeed, using expressions (|l|), (|), (0) and the definition of u(r), we derive 

i 

^ = I- J d^y J $(r) (1 +w(r;7,n))K(r;7,n)c/ 3 r, (24) 



where, according to the notation mentioned above, the identity K(r) = K(r; 7, n) is meant. 
An equation for K(r) can be determined by integrating ([12]) over q, which results in the 
following: 

* 2 K(r) 



2 m 



( V 2 K(r) - (r / . V 2 u(r)) = u(r) + 
V l + u(r) ) 

+ n T (l+u(r)) J $(| r-y |) (1 + u(y)) K(y) d 3 y . (25) 



The boundary conditions additional to (p5|) are of the form 

lim K(r) = 0, K(0) < 00 (26) 

1 — >oo 

and follow from ( P0| ) and (|2"T|). 

Now, let us return to the weak coupling regime to compare our results with the data of 
the Bogoliubov approach. Using (H), (0) and ([TJ|) we arrive at the equality 

\T(q) + ri7$(g)j 
which, after integration, can be written as 

n(q)=f(x)= 1 - flnfl + - * = -^. (27) 



2 V V a;/ 1 + n <&(<?) 

Expression ( |2T| ) should be compared with the result of the Bogoliubov approach that, up to 
the first order, is given p] by 

"(^-^Hvra- 1 )' (28) 

The values of f(x) and fBog{x) for 0.02 < 2 < 8 are listed in Table 1. Besides, for large 
x ^> 1 and small x -C 1 we have 

/(&) ^ /b os (z) - ^ ( x > x ) ' 

/(z) - -^tt > /Bo fl (a:) ^ z-4= (x < 1) . 
2 2v2x 
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So, the model considered provides quite satisfactory estimates of n(q) for x > 0.04 . However, 
it does not yield adequate values of the boson distribution over momenta when x — *■ 1 . This 
feature is not unexpected because the interval of small x corresponds to the region of small 
q . In other words, the behaviour of n(q) in the case q — > is significantly influenced by the 



(see the discussion in [§],[|). The model of Ref.l is able to yield good evaluations for the 
pair correlation function g(r) — 1 when r ^ r c and r ~ r c (here r c stands for the correlation 
length). This means that it can provide reasonable estimates of n(q) for qr c ~ 1 and qr c ^> 1 
and relevant values of the macroscopic quantities such as the mean energy or the condensate 
fraction n /n . Let us, for example, consider the cold many-boson system with the Coulomb 
interaction in the weak coupling regime. Remark that new interest in this system has 
been inspired by curious results of applying the Bose and Bose-Fermi liquid models in 
the investigations of high-temperature superconductors 0. In the case of charged bosons 
$(g) = Ane 2 /q 2 (e is the boson charge) and thus x = 2(q/A) 4 , where A = (IQnnme 2 /Ti 2 ) 1 / 4 
is the inverse screening radius. Using equalities (|I|) and Q27D we obtain 



asymptotics of g(r) — l for r — ► oo which is not reproduced by equation (^j in a proper way 




(29) 



Integration by parts allows one to rewrite expression (|29D in the form 




(30) 



o 



Further, with the equality (see 0) 




oo 



(31) 



one can derive 




(32) 



where r$ stands for the Brueckner parameter given by 




S 



The leading term of the result for the condensate fraction in the Bogoliubov approach |J is 
of the form 



71 



« 1 -0.2114 r| /4 . (33) 

Bog 



The reasonable character of the model considered can also be illustrated with the calculations 
of the ratio of the mean kinetic energy (T) to the mean interaction energy (U) . Using (| 
and x = 2(q/A) 4: , for the quantity (T)/N we have 
(T) V 



N (2tt) 



n(q)T(q) d q 



h 2 A 5 °? / / 1 N 1 N 

Vl + il-T^W (34) 



(4 7 r) 2 2 1 /4 mra 7 V V t 4 / 1 + t 4 
v y o 

Integrating over t by parts and taking account of fl31~|), one is able to find 

i|_4(|)- t ,3/^o, 21 3 t ,-, (35) 



4 

me 



where Ry = ■ Note that when e is equal to the electron charge, we obtain Ry « 13.6 eV . 
The mean interaction energy of the system of charged bosons calculated within the model 



of Ref.l to the first order in the weak coupling approximation, is written as 

-1.107 rj 3/4 . (36) 



(U) _ f3\^ 4 3/4 _ , 1 



NRy \2, 
Equality ( |36"D can readily be derived via the expression 

// \ g 2 If 47re 2 

^(r; 7 , n) - 1 J - d 3 r = lim — J u(q; 7 , n) — d 3 q 



T 

and with (|i~7|) and (|3l|) . Combining (|35D and (^) gives 

ifH =a2 - (37) 

The leading term of the mean kinetic energy in the Bogoliubov approach can be calculated 
with the help of ( |28|) on the analogy of deriving (|34|) 

(r>flos _ »m- 7/1 +** _,W (38) 



AT (4vr) 2 2 1 /4 mn y \^/WTW 



o 



Numerical integration yields 

1 +t 4 







Vt 8 + 2t 4 



1 t 4 dt « 0.2015 



and we hence get 



(T) 



^£ = 0.2008 r 5 3/4 . (39) 



The mean energy of the charged bosons taken up to the first order in the Bogoliubov ap- 
proximation has the form || 

E, 



NRy 

Relations (|39|) and (|40|) lead us to the ratio 

(T) Bo9 <T> 



|2£ « -0.8026 r s 3/4 . (40) 



Bog 



0.2001 (41) 



| (C/) Bog | | -Ebo 9 - (T) B og 

that is in nice agreement with (|37|). 

Thus, considered in the weak coupling approximation, equations (H) and ([12]) provide 
good estimates for the pair correlation function and boson momentum distribution in rele- 
vant regions of distances and momenta. This makes it possible to get adequate evaluations of 
the condensate fraction and main thermodynamic quantities of the cold many-boson system 
with weak interaction between particles. Owing to the correct account of the short-range 
correlations in @ ( see |l]3]), one may expect that equations (H), dl2|) and, hence, (|5|) are 
able to yield reasonable data beyond the weak coupling regime as well. 
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TABLES: 



Table 1 



X 


/(*) 


fBog{x) 


X 


/(*) 


fBog{x) 


0.02 


1.476 


2.037 


1.0 


0.097 


0.077 


0.04 


1.148 


1.320 


1.5 


0.056 


0.046 


0.06 


0.964 


1.008 


2.0 


0.036 


0.030 


0.08 


0.839 


0.824 


4.0 


0.012 


0.010 


0.10 


0.745 


0.700 


6.0 


0.006 


0.005 


0.50 


0.216 


0.171 


8.0 


0.003 


0.003 
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